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We demonstrate that many features ascribed to strong correlation effects in various spectroscopies 
of the cuprates are captured by a calculation of the self-energy incorporating effects of spin and 
charge fluctuations. The self energy is calculated over the full doping range from half filling to the 
overdoped system. In the normal state, the spectral function reveals four subbands: two widely split 
incoherent bands representing the remnant of the two Hubbard bands, and two additional coherent, 
spin- and charge-dressed in-gap bands split by a spin-density wave, which collapses in the overdoped 
regime. The resulting coherent subbands closely resemble our earlier mean-field results. Here we 
present an overview of the combined results of our mean-field calculations and the newer extensions 
into the intermediate coupling regime. 
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Over the past few years evidence has been mounting 
that correlation effects in the cuprates are not as strong 
as previously thought and that these materials may fall 
into an intermediate coupling regime. If so, the transi- 
tion to an insulator could be described in terms of Slater 
physics by invoking a competing phase with long-range 
magnetic order, rather than Mott physics requiring a dis- 
ordered spin liquid phase with no double occupancy. The 
most striking evidence in this direction perhaps is the ex- 
istence of quantum oscillations in underdoped cuprates 
associated with small Fermi surface pockets. [l|, Strictly, 
quantum oscillations have only been observed in strong 
magnetic fields, so that it is possible that the ordered 
phases are field- induced. Q Nevertheless, it is clear that 
a phase with long-range order lies close in energy to the 
ground state. Additional evidence comes from the recent 
calculation of the magnetic phase diagram of the half- 
filled t-t' -U model by Tocchio, et a/.fl They found 
that the ground state is either non-magnetic at small 
U or possesses a long-range magnetic order, except for 
a small pocket of spin-liquid phase at values of U and 
too large to be relevant for the cuprates. Further, 
the high-energy spectral weight associated with the 'up- 
per Hubbard band' decreases with doping too fast to be 
consistent with strong coupling (no double occupancy) 
models. 

There have been a number of recent attempts to ex- 
tend strong coupling calculations into the intermediate 
coupling regime. Yang, et a/.[3] have introduced a phe- 
nomenological self energy that is similar to that of an 
ordered magnetic phase, but with a special sensitivity to 
the magnetic zone boundary. Paramekanti, et a/.[8[ have 
carried out variational resonant-valence bond (RVB)-like 
calculations where the requirement of no double occu- 
pancy has been relaxed. However, given that the ground 
state is close to being magnetically ordered, one could 
also approach intermediate coupling by including effects 
of fluctuations in a weak coupling scheme, where the 
lowest order (Hertree-Fock) solution already describes 
a long-range ordered state. This article presents an 



overview of our ongoing work in this direction. 0, Q 

We begin with a brief recapitulation of the mean-field 
calculations, which by themselves can provide a good de- 
scription of the doping dependence of the low-energy co- 
herent part of the electronic spectrum. The model is 
nearly quantitative for electron-doped cuprates, where 
the competing magnetic order is known to be (w, n) an- 
tiferromagnetic over the full doping range. Such a model 
can describe many features of the hole-doped cuprates 
as well, even though the competing order [s] are known 
to be more complicated. 13, U| Within this model, for 
electron-doped cuprates, the ground state at half-filling 
is an antiferromagnetic insulator, doping simply shifts 
the Fermi level into the upper magnetic band produc- 
ing an electron pocket near (tt,0), and the magneti- 
zation decreases with doping until magnetism collapses 
in a quantum critical point near optimal doping. This 
quantum critical phase transition in fact involves two 
Fermi surface driven topological transitions 12f , the first 
near x = 0.14 — 0.15 where the top of the lower mag- 
netic band crosses the Fermi level producing hole pockets 
near (n/2,n/2), and a second transition near x = 0.18, 
where the hole and electron pockets recombine into a 
single large Fermi surface. The model has been able to 
describe angle-resolved photoemission (ARPES),[l3, 14 1 
resonant inelastic x-ray scattering (RIXS),[15] and scan- 
ning tunelling miscoscopy (STM)[16j spectra, and the 
unusual pairing symmetry transition with doping seen in 
penetration depth measurements Recently, quantum 
oscillations were observed in electron doped cuprates at 
several dopings, showing a crossover from the hole pocket 
at lower dopings to the large FS at the highest doping. [l8[ 
Furthermore, the areas of the FS pockets measured by 
quantum oscillations are well predicted by the model for 
the electron doped case, while hole doped cuprates re- 
main controversial in this aspect. 

Fluctuations modify the above picture in several ways. 
First, in two-dimensional materials, critical fluctuations 
are well-known to eliminate long-range order and drive 
the antiferromagnetic transition temperature to zero in 
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accord with the Mermin- Wagner theorem, so that over 
a wide range of temperatures only a pseudogap remains. 
[The observed Neel order is driven by small deviations 
from isotropic two-dimensional magnetism.] These fluc- 
tuations can be accounted for in a self-consistent renor- 
malization scheme[lj|, and are necessary to describe the 
response of the system at higher temperatures. Fluctua- 
tions also modify the low-temperature physics at higher 
energies, leading to the high-energy kink or the water- 
falls seen in ARPES 0-iiLeffects of the ARPES matrix 
element nothwithstanding (23l427| . We have recently in- 
troduced the quasiparticle-GW (QP-GW) scheme to ac- 
count for these fluctuations. 0, [t| In this way, we have 
been able to describe the doping dependence of the op- 
tical spectra 28-30], including both the 'Slater-like' col- 
lapse of the midinfrarcd peak with doping and the 'Mott- 
like' persistence of a high-energy peak into the overdoped 
regime. The model also quantitatively accounts for the 
anomalous spectral weight transfer to lower energies with 
doping in the cuprates.[6| 

In a GW-scheme, the self-energy is calculated from a 
variant of the lowest-order 'sunset' diagram, a propagator 
dressed by the emission and reabsorption of a bosonic 
operator, 

1 f°° dujp 

S(k, a, iu n ) = - > r) a ,a' / — — r(k, q,w„, w p ) 

2 ^ ■/-«> 2n 
G(k + q, a', ico n + uj p )lm[W aa ' (q, u p )] , (1) 

Here, a is the spin index. W ~ U 2 x denotes the interac- 
tion, which involves the Hubbard U and the su scep tibility 
X in the random phase approximation (RPA)|3l|, and T 
is a vertex correction. [32| -q^y gives the spin degrees of 
freedom, which takes value of 2 for transverse spin and 1 
for longitudinal and charge susceptibility. The model in- 
volves three Green's functions, the bare Go, the dressed G 
given by Dyson's equation G _1 = Gq 1 — S, and an inter- 
nal Green's function Gi n t which will be described further 
below. A number of different variants of the GW scheme 
can be constructed, depending on the specific Green's 
function used in evaluating the G and the W in the con- 



volution integral of Eq. l.[33| Using the bare G in both 



G and W in the so-called 'Go Wo' scheme corresponds to 
lowest order perturbation theory. Using the dressed G 
in both G and W (i.e. the GW scheme) leads to fully 
renormalized propagator and interaction corresponding 
to an infinite resummation of diagrams. However, this is 
still not the exact self-energy because of the missing ver- 
tex corrections. In fact, GW scheme often gives worse 
results than Go Wo version when the vertex corrections 
are omitted. Bearing all this in mind, our approach is an 
intermediate one, in that it is based on the convolution 
of an intermediate coupling Green's function and inter- 
action. In this sprit, we first calculate the self-energy of 
Eq. 1 by using a parameterized G = Gi n t{Z), and then 
calculate Wi n t exactly based on this Gi n t , and determine 



the renormalization parameter Z self-consistently. 

More specifically, we write Go in terms of the unrenor- 
malized LDA dispersion. Note that our tight-binding 
hopping parameters are not free parameters, but are the 
best representation of the first-principles LDA dispersion. 
All renormalizations, giving rise to the experimental re- 
sults, are embedded in the computed S. W is the sum 
of the RPA spin plus charge susceptibilities calculated 
using Gi n t rather than Go- The key lies in the choice of 
Gi n t ■ Our strategy is to construct the best one parameter 
model for Gi n t (Z) with Z chosen to minimize G — Gi n t ■ 
[Gint = G of course yields the full GW.] To motivate 
our choice, we recall that the main effect of £ at low 
energies is to renormalize the dispersion from the LDA 
values to those observed in experiments (e.g., ARPES). 
This renormalization, which amounts to a factor of 2-3, 
is relatively modest in that the mass does not diverge, 
and depends weakly on &;.[34[ That is, approximately 



£k, ARPES = Zck,LDA- 



(2) 



Hence, we choose perhaps the simplest Gint which repro- 



duces this dispersion renormalization, G int 
with £j n j = (1 — Z _1 )cj, so that 

G mt (k, u) 



G n 



u - ZevhDA + iS 



(3) 



The above expression refers to the paramagnetic phase, 
but the extension to a magnetically ordered phase is 
straightforward where G, x> an d S become (2 x 2) tensors 
for the (n, 7r) antiferromagnetic order. [T 
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FIG. 1: DOS computed from Gint is compared with the fully 
dressed DOS with a vertex correction V = 1/Z and without 
the vertex correction (i.e. T = 1) at a representative doping 
of £=0.04 as discussed in the text. 

Self-consistency is obtained by choosing Z such that 
the low energy dispersion is the same for G and Gi n t- 
This is illustrated in Figs. 1 and 2. Fig. 1 compares the 
density-of-states (DOS) associated with Gint (red line) 
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(a) x-0.1 (b) x-0.15 (c) x-0.18 (d) x=0.22 
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FIG. 2: Spectral intensity plots along the high symmetry lines for antiferromagnetic (AF) (x = 0.10), AF+d— wave super- 
conducting (SC) (x — 0.15), d— SC (x — 0.18), and the paramagnetic (PM) (x = 0.22) states. Full self-energy E is included 
in all cases. Antiferromagnetic order parameters are computed self-consistently.[y] Self-energy is computed including the an- 
tiferromagnetic gap, but does not incorporate the low energy superconducting gap. Black lines depict the dispersions that 
enter into Gint, but the corresponding spectral weight of the band is not shown for simplicity 14]. An artificially large value of 
superconducting gap (A = 30 meV) is used so that the effects of superconductivity can be visible on the energy scale of the 
figure. 35] Arrows mark the start of the high-energy kink or the waterfall in the spectrum below the Fermi energy. 



and G, either with (blue line) or without (green) a vertex 
correction in the antiferromagnetic state. As discussed 
further below, G and Gmt— dressed DOSs clearly match 
well with each other in the low energy regime where 
both spectra show the spin density wave dressed upper 
and lower magnetic bands, consistent with our earlier 
mean-field results [3, G3, \$M- At high energies, however, 
Gi n t fails (by construction) to reproduce the incoherent 
hump features associated with precursors to the upper 
and lower Hubbard bands. 

Fig. 2 compares the spectral weights, A — —Im(G)/n, 
with the corresponding dispersions ek,mt (black lines) 
that enter into Gi n t- As noted earlier, the latter provide a 
very good fit to the low energy coherent features for the 
entire doping range in electron doped Nd2_ a; Ce 2 ;Cu04 
(NCCO). Self-consistency ensures that Gi n t provides the 
best approximation to the full G. As an added benefit, 
Z in general renormalizes Gk %nt 

to values close to those 
found in our earlier mean-field studies 14, 17, 35| . The 
obtained self-consistent values of Z decrease almost lin- 
early with doping, as seen in ARPES [36|. In this way, 
the QP-GW model reproduces the results of our earlier 
mean-field calculations in the low energy region, while re- 
vealing new physics at higher energies [e.g., the waterfall 
effect] . 

We now comment on some applications of the present 
QP-GW model. The waterfall physics seen in ARPES 
spectra of the cuprates is a direct consequence of the 
self-energy correction, which introduces a peak in scat- 
tering at intermediate energies below as well as above 
the Fermi level as seen in Fig. 2.[27j This scattering splits 
the spectrum into a low-energy coherent part and a high- 
energy incoherent region. While the near-Fermi-level dis- 
persion changes substantially as the magnetic and the 



superconducting phases evolve with doping, the over- 
all energy regime of the waterfall phenomenon remains 
fairly doping independent (marked by arrows in Fig. 2), 
consistent with experiments. [21| In the pseudogap region 
[x = 0.10, Fig. 2(a)), the resulting 'four band'-like struc- 
ture (two magnetic bands and the two Hubbard bands) 
agrees well with cluster 37( and quantum Monte Carlo 
calculations [3^|. Near optimal doping d— wave supercon- 
ductivity coexists with the antiferromagnetic state in a 
uniform phase 17 , 3|| resulting in further splitting of the 
coherent bands as seen in Fig. 2(b). The coherent bands 
approach the Fermi level with increasing spectral weight 
as the pseudogap collapses at a quantum critical dop- 
ing near x = 0.17 in both the electron and hole doped 
case[l(|. On the other hand, the Hubbard bands move 
towards higher energy as the doping increases and the 
spectral weight associated with these bands decreases, 
consistent with optical spectra. 2^, 13] 

Notably, the lifetime broadening of the quasiparticle 
states originates from magnetic scattering in the QP-GW 
model. This broadening has a non- Fermi liquid form with 
a significant linear-in-T component, particularly near the 
Van Hove singularity (VHS).[1] Furthermore, in describ- 
ing broadening of ARPES features in the superconduct- 
ing state in our earlier mean-field model, we phenomeno- 
logically introduced an elastic small angle scattering con- 
tribution of similar non-Fermi-liquid form[35l l39l . |40j 



E"(w) = sgn(ui)Co 



(4) 



Here Co = 100 meV and ljq = 1.6 eV are determined 
from a fit to the ARPES energy distribution curves 
(EDCs). The exponent p is of physical significance in 
determining quasiparticle character [4l|. We found that 
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FIG. 3: Imaginary part of the self-energy for NCCO (red 
line) is compared with the form used in Eq. [4] Note that 
the model small angle scattering formula was only applied to 
filled states, uj < 0. 



p = 3/2[35, 3{| 4(| applies for electrons as well as hole 
doped cuprates in reproducing the ARPES spectra and 
quasiparticle interference (QPI) pattern seen in scanning 
tunneling spectroscopy. This scattering is particularly 
important in that it allows a finite spectral weight near 
the Fermi level even in the presence of a pseudogap, 
thereby revealing the leading edge superconducting gap 
at all momenta lfl 35| • In Fig- 3 we show that the imag- 
inary part of the self-energy computed from Eq. 1 above 
reproduces this phenomenological form very well in the 
low-energy region, both in magnitude and T-dependence, 
indicating that magnon scattering underlies the anoma- 
lous exponent p = 3/2. 

We emphasize that in our QP-GW model the bare 
dispersion is taken directly from LDA, and the self- 
consistently determined Z renormalizes this into a dis- 
persion which matches the bands seen in ARPES exper- 
iments. Hence, the model reproduces our earlier mean- 
field results, but with fewer parameters, since the dis- 
persion is calculated self-consistently rather than being 
derived from experiment. 

As noted above, the high energy features are absent in 
Gmt ■ This is also clear from Eq. |3] If we integrate the 
spectral weight A int {k,uj) over all frequencies we get Z, 
not 1, so that Gmt accounts for only the coherent part 
of G. The incoherent part of weight, 1 — Z, is thus not 
accounted for. Clearly, this could be done straightfor- 
wardly by including a pair of broadened Lorentzians, but 
this will add additional parameters in the computation 
of the self-energy, not to mention associated vertex cor- 
rections. Therefore, we have chosen to first explore the 
QP-GW model without these complications. This is also 
the reason for calling our approach as QP-GW because 
it focuses on the QP part of the spectrum in evaluating 

Gint- 



In order to better understand the susceptibility Xint, 
a comparison with our earlier mean-field result, xmf 
is instructive. Since Gi n t differs from Gmf only by 
an overall multiplicative factor of Z, Xint differs from 
Xmf by a factor of Z 2 . This has important conse- 
quences. Fermi liquid theory requires that both the dis- 
persion and the spectral weight be renormalized by in- 
teractions, and for a weakly fc-dependent self energy, as 
in the present case, Zd% S p — Z w . This is true for Gmt 
and G, but Z^ = 1 for Gmf, which causes mean-field 
theory to overestimate the tendency towards instabil- 
ity. As the bandwidth decreases (Z —5- 0), the density 
of states and susceptibility must increase as 1/Z, in or- 
der to keep the total electron number fixed since there 
are no incoherent states in mean-field theory. Instability 
is controlled by the Stoner factor, Uxo{q,u = 0) = 1. 
Since Xo(q,0,MF) = Xo(q, 0, LDA)/Z, a small Z en- 
hances the probability of instability. On the other hand, 
for G or Gmt only the low-energy quasiparticle degrees 
of freedom contribute to the instability, as reflected in 
Xo(q, 0, mi) = Zxo(q, 0, LDA). Hence, large fluctua- 
tions (small Z) actually reduce the probability of con- 
densing into any one mode. Equivalently, if we rewrite 
the Stoner factor in terms of the LDA susceptibility, it 
becomes U e f jXo.lda — 1, with f/ e // = ZU . Thus t and 
U should both be renormalized by factors of Z, leaving 
t/U invariant. In our mean-field treatment, we had to 
assume that the effective U was doping-dependent, and 
this Z-correction accounts for part of that doping depen- 
dence. . 

It should be noted that an accurate calculation of the 
susceptibility and the resulting self-energy is fairly com- 
puter intensive as it involves a three-dimensional integral 
(k x , k y ,uj) for xo and a similar three-dimensional integral 
over x's for E. Fortunately, we find that E has only a 
weak fc-dependence, so we need to calculate it only over 
a few fc-points and use the average. Clearly, this is not a 
limitation of the model, and the full A-dependence could 
be calculated. However, this would make accurate self- 
consistent calculations substantially more time intensive. 
We have explored the use of a vertex correction for E, but 
the results are not too sensitive. We have typically taken 
r = 1/Z, which puts somewhat greater weight into the 
incoherent bands as seen by comparing blue and green 
curves at higher energies in Fig. 2. 

The present scheme can straightforwardly incorpo- 
rate the full fc-dependence of the susceptibility based 



on a realistic material specific band structure. This 
is important for delineating the nature of competing 
ordered phases, which are different for electron and 
hole doped cuprates. Moreover, our self-energy pro- 
vides a tangible basis for going beyond the conventional 
LDA-based framework for realistic modeling of various 
highly resolved spectroscopies, providing more discrim- 
inating tests of theoretical models. In addition to the 
ARPES spectra discussed above, a note should also be 
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made in this connection of the optical spectroscopy, [6j 
STM 0, 0, E|L RIXS0, S||, x-ray absortion spec- 
troscopy (XAS) 45|and other inelastic light scattering 
spectroscopies |2a, l46M48j to help piece together a robust 
understanding of the nature of electronic states in the 
cuprates and their evolution with doping. 

In summary, we have shown that our intermediate cou- 
pling model of self-energy, which is based on the spin- 
wave dressing of the quasiparticles, can explain many 
anomalous features of the cuprates. At low energies, the 
model reproduces our mean field results for the coher- 
ent bands in ARPES,[i3] optical, @ and RIXS,[Ii] with 
self-energy corrections renormalizing the large widths of 
the LDA dispersions. At high energies, we obtain the 
waterfall features which represent a splitting off of the 
incoherent bands, precursors of the Mott gaps seen in 
ARPES and optical studies. In the underdoped regime, 
the coherent in-gap bands reproduce both the four-band 
behavior seen in quantum cluster calculations and the 
magnetic gap collapse found in the mean-field calcula- 
tions and a variety of experiments. These results clearly 
suggest that the cuprates can be understood within the 
intermediate coupling regime with an effective U value 
substantially smaller than twice the bandwidth. 
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